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Abstract. Motivated by diffusion processes on metric graphs and ramified spaces, we consider an abstract setting 
for interface problems with coupled dynamic boundary conditions belonging to a quite general class. Beside well- 
poscdness, we discuss positivity, L°°-contractivity and further invariance properties. We show that the parabolic 
problem with dynamic boundary conditions enjoy these properties if and only if so does its counterpart with 
time-independent boundary conditions. Furthermore, we prove continuous dependence of the solution to the 
parabolic problem on the boundary conditions in the considered class. 



Elliptic systems with coupled boundary conditions have been attracting broad attention at least since pQ. A 
classical approach is based on interpreting interface conditions of an elliptic system as boundary conditions of a 
vector- valued elliptic equation. This leads to introducing differential operators acting on spaces of vector- valued 
functions. A parabolic theory for this kind of operators has been recently developed, see e.g. [41 [21]. 

A particularly interesting application of the theory of elliptic systems is given by so-called networks and quan- 
tum graphs, see e.g. [SJ 139) and references therein. Their generalisation to n-dimensional problems has appeared 
already in [43], where the related notion of ramified space has been proposed. Having in mind applications to 
quantum graphs, Kuchment has proposed in |38 a class of coupled, time-independent boundary conditions for 
1-dimensional elliptic systems. Kuchment's formalism allows for a very efficient variational approach, but the 
tradeoff is that his boundary conditions are only a proper subset of those considered in |l a - or, in the specific 
context of quantum graphs, in 37 . However, it is remarkable that Kuchment's conditions give rise exactly to 
all self-adjoint realisations of the Schrddinger operator on a metric graph, under a mild locality assumption. 

In the companion paper |17j . Cardanobile and the author have generalized Kuchment's formalism to the 
case of n-dimensional vector-valued diffusion and characterized several properties of the parabolic problem in 
dependence on the chosen boundary conditions. The aim of this paper is to provide the extension of the theory 
in [17] to the case of dynamic boundary conditions of Wentzell-Robin-type. 

Although we are soon going to consider the general case, let us start by briefly focusing on the 1-dimcnsional 
setting of networks (or quantum graphs). 

Example 1.1. Let N £ N and consider the prototypical case of a diffusion problem 



on a metric graph - more precisely, on a semi-infinite star with N edges ei,...,ejv on whose center a dy- 
namic Kirchhoff-type boundary condition is imposed along with a standard continuity assumption. Each edge is 
parametrized as a (0, oo) -interval, where is identified as the center of the star. Therefore, the function Uj de- 
scribing the diffusion on the edge ej maps [0,oo) x [0,oo) to C, while ip ■ [0,oo) — » C describes the time evolution 
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1. Introduction 



(TDPS) 




u"{t,x), t > 0, x G (0,oo), j = 1, 

u e (t,0) =: ip(t), t>0, j,£=l,...,N, 
EjLi«$(*,0) t>0, 



.,N, 



2 



DELIO MUGNOLO 



of the common boundary value in the center. It is known that the associated initial value problem is well-posed, 
as discussed, e.g. in [3l [lOj [55] . 

Laplace operators with dynamic boundary conditions appear as limiting cases of approximation schemes con- 
sidered in [40[ 126] . The cable model of a dendritical tree proposed by Rail in |62] also leads to analogous network 
diffusion problems, cf. [14U57] : a thorough biomathematical investigation of them has been performed in a series 
of four papers beginning with |44j . 




A semi-infinite star with 6 edges. 

We can rephrase (TDPS) by considering the orthogonal projection Py of C N onto the subspace Y := (1) 
spanned in C N by the vector 

1 :=(!,..., 1). 

Observe that the unknown can be thought of as a function u : (0, oo) — > C , so that the network diffusion problem 
simply becomes 

u(t, x) = u"(t, x), t > 0, x <E (0, oo), 

with suitable boundary conditions in 0. More precisely, the continuity condition in the star's center - given by 
the second equation in (TDPS) - amounts to require that u(t,0) € (1) for all t>0, i.e., 

P Y (u(t,0))=u(t,0), t>0, 

while the dynamic boundary condition equivalently reads 

u(t, 0) = P Y (u(t, 0)) = NP Y (u'(t, 0)) = -Py ( — (t, 0) J , t > 0. 

Hence, the dynamic boundary condition is an equation living in the (1- dimensional) boundary space Y = (1). 

This kind of boundary conditions also arises in the mathematical modelling of string networks with masses at 
the nodes. They play an important role in the control theory of wave and beam equations: investigations in this 
direction go back at least to [41 [ §2.7] and |31j . 

The goal of the present article is to generalize the setting discussed in the above example. Let be a smooth 
open domain in R™ with boundary T :— dfl. Let H be a separable complex Hilbert space. In particular, Bochner 
spaces L 2 (fl; H) and L 2 (T; H) become separable complex Hilbert spaces when endowed with the canonical scalar 
products 

(f\9)mn;H) ■= I (f(x)\g(x)) H dx, f,g e L 2 (n-,H), 
Jn 

and 

- J^f(z)\g(z)) H do-(z), f,ge L 2 (T; H). 
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Let y be a closed subspace of L 2 (T;H) and hence a Hilbert space in its own right with respect to the scalar 
product induced by L 2 (T; H). Vector- valued Sobolev spaces can be introduced recursively just like in the scalar- 
valued case. I.e., one first lets H°(£l; H) := L 2 (fl; H), hence defines 

H k (Q;H) := {/ £ H h -\£l;H) : 3Vf := g e L 2 (n ; H n ) s.t. 
l ' j J n f(x)Wh(x)dx = - J n g(x)h(x)dx for all h e C c °°(fi;C)} , fc = l,2,..., 

and finally introduces spaces of fractional order by standard complex interpolation. (Here we denote by H n the 
Hilbert space defined as the Cartesian product of n copies of H .) In particular, H 1 (Vt\ H) is a Hilbert space with 
respect to the scalar product 

(1.2) U\g)m { n;H) ■■= f (Vf(x)\Vg(x)) Hn dx+ [ (f(x)\g(x)) H dx, f,geH l {£l;H), 

Jn Jn 

We emphasize that vector-valued Sobolev spaces are introduced using scalar-valued test functions, hence the 
integral appearing in (|l.lj) is vector- valued (i.e., a Bochner integral) whereas those appearing in (|1.2|) are scalar- 
valued (i.e., Lebesgue integrals). It is well-known that the usual trace and normal derivative operators 

(7(0) 9 u m. u\ T £ C(T) and C'(fi) 3 « 4 - £ C(r) 

extend to operators acting between Sobolev spaces of scalar-valued functions. In fact, they can be canonically 
defined in the vector-valued case, too - e.g. by means of [30l Thm. 4.5.1]. With an abuse of notation we therefore 
denote by U|p and ^ the trace and normal derivative (in the sense of distributions) of a function u : Q — > H . 

We are now in the position to generalise the one-dimensional setting presented in Example II .11 by allowing for 
more general coupling conditions at the interface and consider the abstract boundary-value problem 

r %u(t) = Au(t), t>0, 
(AS) \ u(t)| r € y, t>0, 

1 »„(t) |r = -Py 9 -^!, t>0, 

where Py denotes the orthogonal projection of L 2 (T; H) onto the closed subspace y. In the 1-dimensional case 
of finite quantum graphs, the investigation of such a problem has been sketched in [341 §4]. 

Example 1.2. Let SI = (0, oo) and H = C N , so that L 2 (T) = L 2 ({0}) = C N . Take 

y:=((l))={ceC N : Cl = ... = c N }. 

Then 

(I ... ^ 

\i ... i, 

and one sees that (AS) is just a reformulation of (TDPS) considered in Examvle \l.l\ 

Example 1.3. Let again = (0, oo) and H = C N . If N = 1 and y — L 2 (T\H) = C, then the first boundary 
condition in (AS) is void and (AS) is the reformulation of a scalar-valued heat equation with Wentzell-Robin 
boundary conditions, see e.g. the recent contributions in [271 1^1 [Ml 165) . If instead y = {0}, then (AS) reduces to 
a heat equation with Dirichlet boundary conditions. For N = 1, these are the only possible choices for y , but for 
N > 2 we have infinitely many new boundary conditions that in some sense interpolate between Dirichlet and 
Wentzell-Robin ones. This is crucial when setting up a Courant-Fischer min-max formula, cf. |llj . 

Example 1.4. For H = <C N the elliptic problem with dynamic interface conditions - a vector-valued version of 
Wentzell-Robin boundary conditions - has been considered in [631 §111.4.5]. In [58| , even more general elliptic 
interface problems have been considerd under the very general assumption that the given system can even consist 
of several metric spaces with different Hausdorff dimensions, see also [12] , see also [3] . 

As already mentioned, the general case of a diffusion equation equipped with coupled (either dynamic or time- 
independent) boundary conditions is mostly motivated by the theories of quantum graphs and parabolic network 
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equations, but it also appears in higher dimensional applications, in particular in biomathematical models - see 
e.g. [48] and references therein. 

In this article we restrict to the case of dynamic boundary conditions only. However, the general case of mixed 
dynamic/time- independent boundary conditions (typically appearing in models from the applied sciences, see 
e.g. [44]) can be easily treated combining the results presented here and those from [17] , 

In Section [5] we introduce our abstract framework and deduce a well-posedness result. The above examples 
suggest that the vector-valued setting - although equivalent to the that based on a network (or ramified space) 
formalism - is more efficient. In fact, its flexibility allows to simply introduce whole families of spaces y. 
Consequently, completely new questions arise. For example, one may wonder how the solution to the heat 
equation with boundary conditions as in (AS) depends on y-. it will be shown in Thcorcm l2.6l that this dependence 
is continuous in norm under very natural assumptions. This result is interesting in that it does not have a scalar- 
valued pendant. We also extend to the vector-valued case a result on continuous dependence on parameters 
obtained in the scalar- valued case in [21] . 

We consider invariance of order intervals and subspaces in Section [3] showing in particular a tight relation 
beween the properties of the heat semigroup governing the problem with time-independent (i.e., Robin- type 
vector- valued) boundary conditions and its dynamic counterpart. We will observe some unexpected phenomena: 
e.g., the semigroups governing these diffusion problems are in general not submarkovian - not even positivity 
preserving. 

To discuss these behaviours in detail, in Section [4] we focus on the setting of Example ll.il It turns out that 
even in the simple context of diffusion on a semi-infinite star with finitely many edges, unexpected dynamics 
arises after chosing appropriate boundary conditions. 

Finally, in Section [3] we briefly discuss the general properties of a similar but different kind of dynamic 
boundary condition, where the normal derivative - rather than the trace - undergoes a time evolution. 

2. Preliminary results 
To begin with, we make our standing assumptions precise. 

As in the previous section, let if be a separable complex Hilbert space, f2 be an open domain in R" with C 1 
boundary T :— dil and y be a closed subspace of L 2 (T; H). In the rest of the paper we are going to investigate 
the general abstract initial-boundary value problem 

M-u{t) = Au{t), t > 0, 



(AV) 



u(t)|r e y, t>0, 
JW% = Py(-^ + (7Ar-<S)n(%), t > 0, 

u(0) = M , 



II 



(0)i r = uq. 



Here 7 G R+, 

SGC(Hi(r;H);L 2 (T;H)) 

and Ar denotes the (dissipative) Laplace-Beltrami operator on the (n— l)-dimensional (differentiable, orientable) 
manifold T, with the convention that 7 = if n = 1, and hence if T only consists of isolated points. The vector- 
valued Sobolev space if 1 (IT; H) can be defined in the usual way as the vector- valued version of the scalar- valued 
space H 1 (r) as introduced, i.e., in [42l §1.7.3]. The Laplace operator appearing in (AV) is defined weakly. 

While weak defining the Laplace operator on open domains is standard, a more detailed introduction of the 
(weakly defined) Laplace-Beltrami operator is in ordeiQ. A definition of the Laplace-Beltrami operator by means 



1 Observe that any differentiable function g : Y — > H is a mapping between the differentiable manifold T and the (trivial) Hilbert 
manifold H, whose tangent bundles are TV = V X R n_1 and TH = H X H , respectively. Accordingly, at any point x £ T the 
derivative Vg(x) : T X T — > T g ( x ^H is a bounded linear operator from R n— 1 to H — hence it can actually be seen as a vector in H" -1 . 
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of Hilbert space techniques has been performed in the recent preprint [8]. In fact, 

(Vr/(-)|Vr0(-))a*-i : T C, f,ge H^H), 
can be defined as the Lebesgue-integrable mapping such that its restriction to any chart (V, £) on T satisfies 

n-l 

where g is the canonical Riemannian metric of the surface V . This expression defines in turn a sesquilinear form, 
and the linear operator associated with this sesquilinear form is the (weakly defined) Laplace-Beltrami operator 
A r . We refer to §1] for details. 

Remark 2.1. Clearly, both A and Ar may be replaced by general elliptic operators with real-valued coefficients 
in pretty much the same way |54j generalizes [9] . Similarly, lower order terms may be added. 

It is known that the right setting for the study of systems of this kind is either the space of continuous 
functions on Q, or else an L p -product space. We are going to follow the latter approach throughout this note. 

Lemma 2.2. The space 

(2.1) Vy := jf := (J) GH^Q-H) x (H s (T; H ) n 3>) j 

is dense in C y := L 2 (Q; H) x y for all s > 0. 

In no confusion is possible, in the following we will write C 2 instead of Ly. 

Proof. This is a slight modification of [541 Lemma 5.6]. More precisely, the assumptions in [54[ Lemma 5.6] 
can be weakened by merely assuming that i? 1 (r; H) n 3^ is dense in the range of the trace operator, instead of 
coinciding with it. This density condition is satisfied by assumption, hence the claim follows. □ 

In the following we set either s = 1 if 7 > 0, or s — \ if 7 = 0. Accordingly, 



Vy 

or 

Vy: 



f = eH\n ; H) x (H\r ; H)ny)^ if 7 >0 
etf^Ojff) x (#3(r ; #)n;yH if 7 = 0. 



/ 

/|r 

We consider a form (as,Vy) defined by 

a s (f , g) := f (Vf(x)\Vg(x)) Hn dx + 7 / (V r /(z)|V r g(z)) H „- 1 ^(z) + (5/|r|fl|r)y. f,g6 ^y, 

where the second addend on the right hand side corresponds to the Laplace-Beltrami operator on the Riemannian 
manifold T (recall that by convention 7 = whenever n — 1). We remark that 

(<5/|r|ff|r)y = (Sf\r\Py9\r)y = (PySf\ r \g\r)y for all f,g e Vy, 

so that the third addend in the definition of (ag, Vy) is well-defined. 

By a principle presented in 17, Appendix] and based on [30., Thm. 4.5.1], the classical Maz'ya inequality 
(cf. §4.11.2]) can be extended to the vector-valued case. Accordingly, in either case is a Hilbert space 
with respect to the norm defined by 



(f \g) Vy := / (Vf(x)\Vg(x)) H ndx + / (V r f(z)\V r g(z)) H n-i da(z) + (f\ r \g\r)y if 7 > 



or 



(i\s)v y := / (Vf(x)\Vg(x)) H ndx + (f\ T \g\ T )y if 7 = 0. 
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Theorem 2.3. The operator Ay t s associated with (as, Vy) generates an analytic semigroup (e tAy - s )t>o with 
angle \ on C 2 . 

The operator Ay^s is dissipative if the operator 5 is accretive and in this case the semigroup (e tAy - s ) t >o is 
contractive. The operator Ay^s is self-adjoint if and only if the operator 5 is self-adjoint and in this case the 
semigroup (e tAy - s )t>o is self-adjoint. The operator Ay.s has compact resolvent if and only if Q,T have finite 
measure, provided that H is finite dimensional: in this case the semigroup (e tAy < s )t>o is compact. 

The proof is based on the approach presented, e.g., in 23, Chapt. VI]. We borrow our terminology from [BJ. 

Proof. We are going to show that (as,Vy) is associated with an operator that generates a cosine family with 
phase space Vy x C 2 in the sense of [3 §3.14]. To this aim, we show that for all 7 G M + the densely defined 
sesquilinear form (as, Vy) is continuous and elliptic (with respect to C 2 ), i.e., 

Rea s (f, f) + w||f ||| a > a\\{\\ 2 Vy for all f G Vy 

for some a > and a suitable wei 

Continuity follows from the Cauchy-Schwarz inequality. Ellipticity (with respect to L 2 ) follows from ellipticity 
(with respect to L 2 (£l; H) and L 2 (T; H)) of the forms associated with the Laplace and Laplace-Beltrami opera- 
tors, corresponding to the first two addends of (as, Vy). The third addend in the definition of as is sesquilinear 
and defined on H 2 (T; H) xH 2 (T; H), hence it can be neglected by a perturbation argument (see [HJ Lemma 2.1]). 
Finally, because 

|Ima s (f,f)| = M5/| r |/| r ) y |< 1151111/11^1 ^H/irlU^-H) 

< -^l|5||||/||ffi(fi ; _f/)||/|rl!L 2 (r;,ff)7 

for some M > and all f G Vy due to boundedness of the trace operator from i7 1 (i7;i?) to H^(T;H), the 
announced generation of a cosine family follows by [22j Thm. 4] . It is known that generators of cosine operator 
functions also generate analytic semigroups with angle j, see [3 Thm. 3.14.17]. 

Because the forms associated with the Laplace and Laplace-Beltrami operators are accretive, accretivity of 
(as, Vy) is clear provided S is accretive. A direct computation shows that (as, Vy) is symmetric if and only if S 
is self-adjoint. The assertion on compactness follows from the Aubin-Lions Lemma, see |641 Prop. III. 1.3]. □ 

The proof of the following is based on (9j Rem. 2.2]. 

Proposition 2.4. Assume to have C 2 -boundary. For all 7 G M+ and S G C(L 2 (T\ H)) the operator Ay t s 
associated with (05,^) is given by 

D(A y>s ) = jf := G Vy : Af G L 2 (n ; H), A r /|r G L 2 (T:H), and |£ G L 2 (T;H) 

A 

d 
>dV 

hence (e tAy - s ) t >o yields the solution to (AV). If in particular f G D(A y , s ), then f G Hi(Q; H) H Hf oc (VL;H). 

Observe that in general Ay would not operate on C 2 if we would drop the term Py. 
Proof. By definition, the operator associated with (as, Vy) is given by 

D(By,s) := {f G : 3g G £ 2 s.t. a (f, h) = (g|h) £2 VheVy}, 
By^sf ■= -g- 

By the perturbation theorem of Desch-Schappacher (cf. [25]) a relatively bounded perturbation does not affect 
the domain of an operator. Hence we can assume w.l.o.g. that 

5 = 0. 



y ' S \- p y£ Py(jA r -S) 
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In order to prove that Ay t s C By, 5 take f , h £ Vy. By the Gaufi-Green formulae and the (weak) definition of 
the Laplace-Beltrami operator we obtain 

fls(f,h)= f (Vf(x)\Vh(x)) H ndx + ~/ f (Vvf\r^Wrh\ r (z)) da(z) 
Jn Jr 



(Af(x)\h(x)) H dx 

^ V ^/(*)IV(*)) ^)-7j^(Ar/|r(«)l5|r(«)) H M«) 
= - J^Af(x)\h(x)) H dx + (yLf( z )\h(z)j - 7 (A r /,r|5|r)3, 

and the operator Ay s has the claimed form. 

Conversely, let f £ D(By j s). The above computation also shows that A/ and Ar/|r are well defined elements 
of L 2 (Q;H) and L 2 (T;H), respectively, and that / has a weak normal derivative in L 2 (T;H). We deduce that 
/ £ H%(Q; H) by [42j Thm. 2.7.4] - suitably extended to the vector-valued case by virtue of [30l Thm. 4.5.1]. □ 

Remark 2.5. The vectors in D(Ay) also satisfy the additional boundary condition 

(2.2) (Au) lr £y and (Ait)|r + Py — + Py {Su\ T - 7A r W|r) = 

for all z £ V. Conditions can be interpreted as a formulation of Wentzell-Robin boundary conditions which 
is stronger than the dynamic one that is usual in the context of L p - spaces. Due to the regularising effect of the 
analytic semigroup (e tAy - s )t>o, these additional conditions are satisfied by the solution (AV) for any time t > 0. 

Consider a sequence (y n )neN 01 closed subspaces of L 2 (T;H) such that the associated sequence of orthog- 
onal projections (Py n ) n eN converges in operator norm. Then its limit is also necessarily a projection and a 
contraction, i.e., an orthogonal projection - say, onto a subspace y. Consider moreover a sequence (iSn)neN in 
C(H^ (T; H); L 2 (T; H)) that converges in operator norm to some S £ C(H^ (T; H); L 2 (T; H)). Now, it is quite 
natural to conjecture that Ay nt s„ converges to Ay t s in a suitable sense. 

Observe that no kind of convergence from above or below of the form family (as n , Vy n ) n eN holds - in our case 
one typically has Vy n n Vy m = Vy for some lower-dimensional J^o, whenever n 7^ m - so that in general Vy is 
not dense in any Vy n . Furthermore, the operators Ay ni $ n and Ay,s act on L 2 (fl;H) x y n and L 2 (£l;H) x y, 
respectively, i.e., they generally act on different spaces. All in all, it seems that well-known results for convergence 
of operators associated with forms (e.g., those due to Kato and Simon) cannot be applied to our setting. Some 
results on approximation of operators acting on different spaces have been recently obtained by Ito and Kappel 
(see e.g. [351 Chapt. 4]), but it seems that they fall short of fitting our framework, too. 

The different approach proposed by Post in |60) and further developed in |53j seems to be more appropriate. 
In order to apply Post's results, we need to impose a structural assumption on y that will prove a significant 
simplification in our framework. 

Theorem 2.6. Let (Y n ) ne jq be a sequence of closed subspaces of H. Consider a further closed subspace Y of H 
and a family (J^ n ) n ef-i °f unitary operators on H such that J^ n Y n — Y for all n £ N. Assume furthermore that 
limn^oo J^ n — Id in operator norm and consider the spaces 

(2.3) y := {/ £ L 2 (T; H) : f(z) £ Y for a.e. z £ T}, and 

(2.4) y n := {/ G L 2 {T; H) : f{z) £ Y n for a.e. z £ T}, n £ N. 
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Let (S n )neN be a sequence of accretive bounded linear operators on H that converges in operator norm to some 
S G £(H) and define linear operators S n ,S € C(H^ (T; H)) by 

S n g :— S n o g, n <E N, and Sg := S o g, g£H^{Y\H). 

Then both families (i?(A, A^s^)^^ and (e*^™' 5 ™ )„ e N of bounded linear operators on Cy converge in operator 
norm to the bounded linear operators R(X, Ay,s) o,nd to e tAy - s on Cy, for all ReA > and for all t > 
respectively. Moreover, if H is finite dimensional and Cl,T have finite measure, then the (discrete) spectrum of 
Ay n ^s n converges to the (discrete) spectrum ofAy $. 



Remark 2.7. Observe that the phenomenon observed in Theorem \2.6\ is intrinsically related to the vector-valued 
case. If in fact dim H = 1, then each sequence (Yn)neN of subspaces of H = C such that (-Py„) n eN converges is 
eventually constant - with value either {0} or H - so that the assertion becomes trivial. 

The proof is based on an abstract convergence scheme discussed in |60[ Appendix] , which we briefly recall for 
the sake of self-containedness. The following collects results from [BD1 Thms. A. 5 and A. 10] 

Proposition 2.8. Let 'H,'Hi,'H,'Hi be Hilbert spaces such that T~i\ <—t % and T~i\ >• % with dense embeddings. 
Let f) : "Hi x "Hi — > C and f) : T-L\ x Hi — > C be continuous, accretive and elliptic (with respect to % and 
H, respectively) with associated operators 21 and 21. Consider operators J € C(H,'K), J € C(H,1-l), J\ £ 
C(Hi,'Hi), J\ £ CfflijHx). Let moreover the above spaces and operators satisfy the following conditions: 



(2.5) 

(2.6) 

(2.7) 

(2.8) 

(2.9) 

(2.10) 

(2.11) 

(2.12) 

for some S > 0. Then 
for some M > 0. 



\\Jf-Jii\\n < S 

WJn-JmWu < S 

|G7f|u)rt-(f|.7u)„| < 5 

- f)(f|Jiu)| < 5 

\\f-JJf\\n < 5 

\\u-jjuWz < S 

WJfW* < 2 

WJHn < 2 



u 



Mi' 



f||M||u|| K > 

flkilluii^^ 

f||Mi, 



u 



Mi' 

Im, 



u 



|i?(A,2l) - JR(\,%)J\\ < MS 



We emphasize that the convergence assertion is rather poor at a numerical level, but fairly strong at a 
functional analytical level: it states convergence in operator norm, rather than just strong convergence as done 
e.g. by the various Trotter-Kato-type theorems. We are now in the position to prove Theorem [ 



Proof of Theorem 



Fix n £ N. We apply Proposition 

•2 

■y 



1~L '. — jOy , 1~L 



| setting 
Vy n , Hi :=V y , 
along with 

f) := (a Sn ,Vy n ) and t) := {a s , Vy). 
Observe that accretivity of f), f) follows from accretivity of the operators S n , S. Define moreover J £ C(H, H) by 



(2.13) 

and 

(2.14) 



Jf:= 



J in o fl 
J in o h 



e H, 



O Ul 

o u 2 



u := 



e H, 
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and moreover Ji := J and J\ := J. 

It is apparent that (|2.5|) . (|2.6p . ()2.11|) and ()2.12j) are trivially satisfied for 5 large enough (and getting smaller 
and smaller as n increases). Moreover, (|2.7I) . (|2.9[) and (|2.10[) hold because J, J are unitary with J* = J. 

Finally, observe that by (|2.3p . (|2.4p . (|2.13[) . (|2.14[) the operators J\, J\ do not depend on space, hence they 
commute with the local operators associated with the forms f), f). Furthermore, for all f € H\ and all u e "Hi 

(5^/| r |«| r )y - (S n f\r\(J in rWhn = {(SJ in ~ J in S n )f lr \u lr )y n 

which converges to because J^ n _1 <S„ J^™ converges to S in operator norm. We conclude that (|2.8[) is satisfied. 

Then the convergence of (R(X,Ay n )) n ^ follows from Proposition 12.81 The remaining assertions follow 
from [501 Thms. A. 10 and A. 11]. □ 

Remark 2.9. Let us consider the case of a more general diffusion equation of the form 

§- t u(t) = V-(Z>Vu(f)), t>0, 
«(t) ]r € y, t>0, 
(AVd) { &u(t)| r = -f^*£& + ( 7 A r -S))u(f)| r , t>0, 

Tt(0) = U , 

u(0)| r = u , 

where D € C 1 (fi; C(H n )) satisfies for some /i > t/ie following ellipticity condition: 

Re(D(x)£\t) H » > MII£IIh» /or a// £ € ff" and aZ/ ieSl. 



XTie subspace y as well as the operator S are now fixed. Then, a variational approach can still be pursued, after 
introducing suitable weighted Bochner spaces Ly D as it has been done in |54j . Due to uniform ellipticity, the 
coefficients do not degenerate on the boundary, yielding that \\-\\c^ D an d II ' II are equivalent norms on Cy D . 

Now, consider a uniformly elliptic family (Dk)keN C C (Ct;jC(H n )) of coefficients such that Dk(x) is self- 
adjoint for all x £ f2 and all fegN. Consider the sesquilinear form ak arising from the problem (AVu fc ) ; k € N, 
whose domains all coincide with Vy . Denote by Afc the associated operator. These operators are uniformly 
sectorial - actually, all their numerical ranges are contained in the negative halfline. If the sequence {Dk)keN 
converges strongly, then (<ifc(f , f))feeN * s a Cauchy sequence for all f G Vy. Therefore, by a known result due to 
Kato (see [35, §VIII.3] y ), (R(X, Ay k ))k£Vi converges strongly for all ReA > 0. By simple functional calculus argu- 
ments this also implies strong convergence of (e zA,E ) fegN j or a n z j n £/j g open right halfplane. This is comparable 
with [12 Thm. 3.1]. A similar assertion concerning convergence in operator norm can also be obtained applying 
Proposition [ 



3. LATTICE-BASED INVARIANCE PROPERTIES 

This section is devoted to the characterisation of qualitative properties of (e tAy ' s )t>o- These can often be 
discussed in terms of invariance of relevant subsets of the state space C 2 - most notably, order interval^. By 
Ouhabaz's well-known invariance criterion, such invariance properties can be characterized by simple, almost 
linear algebraic properties of a quadratic form. In a more general form presented in |45l Thm. 2.1], Ouhabaz's 
criterion can be stated as follows. 

Lemma 3.1. LefH be a separable Hilbert space and a. a sesquilinear form with dense domain V that is continuous 
and elliptic with respect to %. A closed convex set C of % is invariant under the semigroup associated with a if 
and only if V is invariant under the orthogonal projection V of H onto C and moreover Re bl(Vu, u — Vu) > 
for all u G V. 



2 It has been observed in 1 191 §5] that also invariances of some subspaces of the state space often reveal important properties of 
the evolution equation. In fact, all results in this section also apply when order intervals are replaced by subspace — of course, even 
dropping all lattice assumptions. 
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In the remainder of this section assume for simplicity that 7 > 0, i.e., 

Vy = {f€ H\n-H) x H\T-H) : /, r G y) . 
(Still, all assertions hold true in the case 7 = with obvious, minor modifications in the proofs). 

To warm up, we start by characterising reality of (e tAy s )t>o. A function in C 2 is called ifR-valued if it takes 
values in the real Hilbert space i2jj underlying H for a.e. x G f2 © T. As a direct consequence of locality the 
forms associated with the Laplace and Laplace-Beltrami operators we obtain the following. 

Proposition 3.2. The semigroup (e y ' s )t>o leaves invariant the real part of C 2 , i.e., the set of all H^-valued 
functions in L 2 , if and only if (SYieg\hiig)y £ K for all g £ y , hence if and only if 

S{f G H*(T;H) : f{x) G H R for a.e. x G T} C {/ G L 2 {T;H) : f(x) G H R for a.e. x G T}. 
In typical applications the space H is a Hilbert lattic^l - hence we will assume henceforth that 

H = L 2 {Z;C) 

for a suitable finite measure space 3, cf. [47l Cor. 2.7.5]. In particular, the scalar products of L 2 (il;H) and 
L 2 {T;H) read now 



{f\g)L 2 (n-,H) : = J J f( x i£,)9{x,€)d£dx and (f\g)L 2 (r-,H) ■= J J f(x,^)g(x,^)d^dx, 

respectively. We can define the positive and negative parts and the absolute value of functions in L 2 (fl;H) 
pointwise, exploiting the lattice structure of H: if u G L 2 (il; H) = L 2 (Q, x 5; C), then u + ,u~ are the function 

Sl3i4!i+(i,.)eL 2 (H;C) and Sl3wH!i»eL 2 (H;C), 

respectively, where u + (x, ■), u~(x, •) are well-defined elements of H = L 2 (S,C) because u(x.-) G L 2 (S;C). 
Observe that the orthogonal projection "P + onto the positive cone of L 2 (£l; H) acts on any u £ L 2 (fl; H) as the 
composition P + o u, where P+ is the orthogonal projection onto the positive cone of H. 
Let a, b £ L 2 (fl; H) = L 2 (Vl x 5; C) and consider the unbounded order intervals 

[a,+co) L 2 (n . H) := {/ G L 2 {fl:H) : a{x) < f{x) for a.e. x £ Q.) 

S {/£L 2 (!!x 5; C) : a(x, < f(x, for a.e. (x, ()el]xS}, 

(-00, 6]i2( n;H ) := {/ G L 2 (fi; ff) : f{x) < b(x) for a.e. x £ il} 

{/£i 2 (!lx 5; C) : /(x, < &(a;, for a.e. (x, £) G Q X 3} . 

These subsets of L 2 (£l;H) are closed and convex. Similarly, for c,d £ L 2 (T;H) = L 2 (£l x 3; C) one considers 
the unbounded order intervals [c, +00)3;, (— 00, d]y. 

Lemma 3.3. Let u, v £ H 1 (CI; H) and u, v G H 1 (T; H) . Then also max{w, v} £ H 1 (fl; H) as well as max{u, v} £ 
i? 1 (r;iJ). Furthermore, 

(3- 1 ) Vmaxju, «}(-,£) = l {u( ., 5) >„ ( .^ )} Vu(-, ^ + l{u(-,£)<<,0} Vw G; 0, 

(3-2) V(t* -»)-(•,£) = l w ., o< „ ( ., e)} V 1 ;(-,0 5 

(3-3) Vmaxju, £}(•,£) = 1 { , 5( . :5) > C( .^ )} Vu(-, + l{«(-,c)<ii(-^)} Vu(-, C), 

(3-4) V(u-v)-(-,0 = 1{uU)<v(-,0}^H;0, 

for a.e. £ £ 3. 

Observe that (l3.1l) . ([3T2)) . (l3.3|) . (|3.4p represent equalities of functions in L 2 (f2;C) and P 2 (r;C), respectively. In 
particular, each "slice" it(-,£) defines a scalar- valued function on f2: it is the differential of this slice-function 
that is denoted by Vu(-, £). The same is valid for u, u, v. 



^ For the necessary notions from the theory of Banach lattices we refer to |56l 1471 . Consequently, also L 2 (Cl;H), y and C 2 are 
Hilbert lattices. Whenever we refer to an operator on a Hilbert lattice as "positive", we always mean "positivity preserving". 
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Proof. The proof goes in several steps. We will repeatedly use the fact that 

u e L 2 (fl; H) ^ L 2 (0; <C) ® H = L 2 (fL; C) <g) L 2 (E; C) = L 2 (n x 5; C) 

in order to reduce a vector-valued relation to a collection of scalar- valued ones: this follows from the elementary 
theory of Hilbert tensor products. 

1) First of all, we recall the following vector- valued extension of [131 Prop. IX. 3], observed in jTTJ Appendix 
A]: Let G : H — >■ H be a Lipschitz continuous mapping and f £ if 1 (57; H). If G(0) = 0, then Go / G H 1 ^; H). 
The proof is an easy modification of [131 Prop. IX. 3]. In particular, this result applies to the case where G is an 
orthogonal projection onto an order interval 

- [a, +oo) L 2 (n . ff) := {/ £ L 2 (57; H) : a < f(x) for a.e. 

(-oo,0 £a(n .H) := {/eL 2 (Q;F):/(z)</?fora.e. 

for a, (3 £ H with —a, (3 £ H + , so that these order intervals actually contain 0. In particular, if / £ 
thenf+J-£H 1 (n-H). 

2) Observe that although / + is formally given by the composition of a Lipschitz continuous mapping on H 
and a function in H 1 ^; iJ), providing a chain rule is not trivial as Rademacher's theorem fails to hold in infinite 
dimensional spaces and it is in particular not easy to understand in which sense the orthogonal projection of H 
onto H+ is "differentiable a.e." , as one would expect in the finite dimensional case. 

To this aim, let / £ H 1 ^; H). By 1), one has in particular and by definition of H 1 ^; H) that 

/ f + (x)Vh(x)dx = - [ Vf + (x)h(x)dx for all h £ C c °°(51; C) 
Jn Jn 

in the sense of if "-valued Bochner integrals. In other words, the above integrals define an element of L 2 (S, C) n . 
Accordingly, 

(J f+{x)^h{x)dx^j (£) = - Q V 7+ (x)h(x)dx) (£) for all h £ C c °°(fi; C) and a.e. £ e S, 
and therefore 

/ f + {x,^)Vh{x)dx = - [ Vf+(x,£)h(x)dx for all h £ C£°(fi;C) and a.e. £ e S : 
Jn Jn 

this can be checked by first considering step functions and then going to the limit. We deduce that f + (■,(,) £ 
iJ 1 (f2;C) for a.e. £ £ E. Since this is a scalar- valued function, we can apply the usual differentiation formula 
and deduce from [29l Lemma 7.6] that 

(3.5) V/+0, - l { /(. ) o>o}V/(-, for a.e. £ £ E. 

Now, because / + £ if 1 (57; H), the weak derivative of / + is necessarily given by (|3.5j) outside a subset of 57 x S 
of zero measure. 

We emphasize that the characteristic function is defined by means of subsets of 57 such that some inequality 
is satisfied by a sca/ar-valued function. In fact the two subsets {/(•,£) > 0}, {/(•, £) < 0} define a partition of 57 
for a.e. £ £ S. 

3) We are now in the position to prove the main assertion. Since u — v £ H 1 (57; H), we deduce from 2) that 
(u - «)+(•,£)> ( u ~ G -ff x (n;C) and the identities 

v(«-«)+(-,o = i{«(.,o>«(.,0}(v«- v«)(-,o. v(« -«)-(-, o = i {u( . >0<<)0} (Vu-vv)(-,o 

hold for a.e. £ € S. Accordingly, both 

-P(-oo,d] u = min{ii, u} = it — (it — w) + and Pf-u.+oo)' 1 ' = max{u, v} = v + (u — v) + 

belong to F 1 and ([3T|) follows. 

The remaining assertions are proven likewise. □ 
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Theorem 3.4. Let a G x S;C) be such that a = (a, air) G Vy and consider the unbounded order interval 

[a, +oo) £2 := [a,oc) L 2(Q. H ) x [a\ r ,oo) L 2^ r . H) 

= {/ G L 2 (n x S; C) : a(x, £) < f(x, £) /or a.e. (x, G x S} 

x {a G L\T x 3; C) : a(z, < for a.e. (z, ^TxS} 

TTien (e tAy - s )t>o leaves invariant [a, +00)^2 if and only if 

(i) Py[a., +oo)£2 C [a, +oo) C 2 and additionally 

(ii) the inequality 



> Va(x,£)(Vf-Va)(x,£,)dxd£ 

Js J {a(-,£)>f (-,$)} 



+7 / / Va{z,0(Vf(z,0 -Va(z,0)da(z) c%+ {Smax{a lr J\r}\(f\r - a lr )-) y 

holds for all f G i? 1 (f2 X S; M) such that f = (/, /| r ) G Vy. 

Proof. By Lemma 13 - 1 1 (e tAy ' s )t>o leaves invariant the order interval [a, oo)£2 if and only if the associated 
orthogonal projection P[ a ,c») 2 l eaves invariant Vy and moreover a(P[ a +oc ) £2 f , f — P[ a ,+oo) c2 ^) — f° r au ^K - 
valued f G Vy. By Lemma T3.31 the first condition is satisfied if and only if P[a.,+oc) y y C 3^. By [351 Lemma 2.3] 
this is equivalent to Py [a, +00 )y C [a, +oo)y. 

The second criterion can be deduced applying Lemma 13.31 and observing that for all f £ Vy 



as(P[*,+oc) c2 f,f ~ -P[a,+oo) £2 f) = -as(max{a,f},(f-a) ) 

= - / (l{a>/}Vo + l {o < /} V/|l {o >/}(V/-Vo)) HB da! 
Jn 

~yj ( 1 { Q |r>/ l r} v ra|r|l{a |r >/ |r }(V r /|r - V r a|r)) ffn _ 1 <M*0 



7 i ( 1 {a |r </ |r }Vr/|r|l{a |r >/ |r }(V r /|r - V r a|r)J <M*) 



max{a|r,/|r}|(/|r 



-°|r- 



(l{a(.,e)>/(-,O} Va ( :E '0 + 1 {a(-,?)</(-«)} V /( a; '0) 



r 



1 {a(-,«)>/(-.«)}( V / - Va )0;0)d£ ^ 
1 {a|r(-,«)>/|r(-,?)} Va |r( z ^), 



, ( 1 {«|r(-,e)>/|r(-,O}( V r/(*>0 " V r a(z,£))d£ M*) 
" 7 // ( 1 {« l r(-,?)</ l r(-,C)} V r/(2 I 0) 



v 1 {a|r(-,0>/|r(-,C)}( V |r/( 2; ; ~ V| r a(2,0)rf^o-(z) 
(<Smax{a| r ,/|r}|(/| r -air)")^. 
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Applying Fubini's theorem we obtain 

a S (P [sL , + oc) c J,f-P[ a .+oo )c J) = - [ [ Vo(s,0(V7- Va)(x,£)dxd£ 

-iff Vo(«,0(V/(«,0 - Va(z,£))d(T(z) dt; 

JsJ{a ir (-£)>f\A;()} 

- (5max{a| r ,/| r }|(/| r - a\ r )~) y . 
This concludes the proof. □ 

Analogous assertions hold for the order intervals (—00, b]^2. 

In general, condition (ii) in Theorem 13.41 will rarely be satisfied. An easy, yet relevant special case is clearly 
that of constant a, i.e., a(x, £) = a for some a £ H and a.e. (x, £) £ f2 x H. In this case, condition (ii) reduces 
to the condition 

(3.6) (5max{o| r> /|r}|(/|r-a|r)~) y <0. 

Observe that if in addition S £ £(L 2 (T; H)), then by Lemma l3~T1 the validity of condition (ii) in Theorem [33] is 
equivalent to the invariance of [a, +oo) L 2( r . ff ) under the semigroup generated by S. E.g., positivity of the semi- 
group corresponds to invariance of [0, oo) C 2 , while £°°-contractivitjfl can be formulated in terms of simultaneous 
invariance of both order intervals [— 1,00)^2, (—00, l]c 2 - 

For the sake of further reference we introduce the following locality assumptions. 

Assumptions 3.5. There exist a closed subspace Y of H , a closed convex subset Ch of H and an operator 
S £ C(H) such that 

• y = {f £ L 2 (T;H) : f(z) £Y for a.e. z£T}. 

• C L 2 {n . H) := {/ g L 2 (n;H) : f(x) £ C H for a.e. x £ 0}, 

• Cy.= {f£y: f(z) £ C H for a.e. z £ T}, 

• C C 2 := C L 2 (n . H) x Cy and 

• Sg = Sogfor all g £ H?(T;H). 

Moreover, £ C or else both 17 and T have finite measure. 

Observe that under Assumptions 13.51 the abstract problem ( AV) becomes a parabolic problem with dynamic 
boundary condition 

(>0,ie!l, 
(>0,zer, 
S)u(t,z)) , t>0, z£T, 

X £ fl, 

z£T. 

If f2 = (0, 00) x R Tl_1 , it is common in the literature to refer to this problem as "diffusion on an open book" (with 
dynamic boundary conditions). If n = 1, this is nothing but the semi- infinite star considered in Example ll.il 

Under the Assumptions 13.51 y is a closed subspace of L 2 (T; H) and C L 2(Q.m,Cy,C£2 are closed and convex 
subsets of L 2 (£l; H), y, and £ 2 , respectively. With an abuse of notation we then write y = Y, S = S and A.y,s 
instead of A.y,s- 



By this we mean contractivity with respect to the norm of L°°(r2 xE;C)x L°°(T X H;C). 



(DBC) 



Q^Uyt, X 

u(t, z 
§- t u{t,z 
u(0, x 
u(0, z 



— Au(t,x), 

e Y, 

= Py(-^(t,z) + (jA r 

= uo(x), 

= V (z), 
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It is crucial that whenever Assumptions 13.51 hold the orthogonal projections of £ 2 (0; H) onto Pl 2 (q-,h), of y 
onto Py and hence of C 2 onto C^ satisfy 

P c L 2 (n , H) f = Pc H of for all / e L 2 (tt;i/), 

p c y 9 = Pc H ° .9 for a11 9 G ^, 



The fact that the projections onto the above subsets of vector-valued function spaces are the compositions of 
a Lipschitz continuous mapping (namely, the projection Pc H ) and a function of class H 1 permits to apply the 
version of a chain rule obtained in Lemma 13.31 Furthermore, due to the local structure of the sets C^q.^ and 
Cy, one sees that in particular 

( P c L 2 (n:H) ° /)|r = {Pc y ° /|r) for all f e Vy. 

Theorem 3.6. Impose the Assumvtions [375l Then C^i is left invariant under (e tAy - s )t>o if and only if 
(i) the inclusion PyCn C Cu holds and additionally 
(ii) the semigroup {e~ ts )t>o leaves Ch invariant. 

Comparable results have been obtained in the context of networks in [36j [34] . 

Proof. First of all, we show that the inclusion Pc c2 Vy C Vy holds if and only if the inclusion PyCy C Cy holds. 
Orthogonal projections onto closed convex subsets of a Hilbert space are Lipschitz continuous mappings, hence as 
already observed by [T71 Lemma 7.3] Pc c2 maps H 1 ^!; H) x H 1 (Y; H) into itself- i.e., the weak differentiability 
conditions is satisfied independently of the boundary conditions. Consequently, Pc 2 Vy C Vy if and only if 
/|r G y implies Pc y f\r G y, for all f € H 1 ^; H). The proof can be completed reasoning as in [XTl Prop. 4.2]. 
By Lemma [33] invariancc of under (e tA ^ s ) t>o is now equivalent to PyCy C Cy and 

Rea s (P Cc2 f , (/ - P Cc2 )f ) > for all f e Vy. 

Due to locality of the forms associated with the Laplacian on f2 and the Laplace-Beltrami operator on Y (and 
hence of the form (as, Vy)), a direct computation shows that 

Rea s (P Cc2 {, (I - Pc c2 ){) = Rc(S P Cy f\ r \(I - Pc y )f\r)y- 

By density, the latter term is > for all f € Vy if and only if 

Re(SP Cy g\(I - P Cy )g)y > for all g e y. 

By a localisation argument this is equivalent to asking that 

Re(SP CH x\(I - P Ch ))x) h > for all x € H. 

A further application of Lemma 13 . 1 1 concludes the proof, since (S ■ \ -)h is the form associated with —5. □ 

In the previous theorem, it is not too restrictive to consider sets of the form Cl^^h) x C^rp.^ - i.e., to 
restrict ourselves to study invariance of sets of those functions pointwise belonging to the same subset of H , both 
on and on the boundary Y. In fact, the following holds. 

Proposition 3.7. Let C, D C H be closed convex subsets. IfCjjir^jjs x D L2 ^. H ^ is invariant under (e tAy - s )t>o, 
then C = D. 

Proof. We only consider the case of fi, Y with bounded measure. The general case will then follow by localisation 
arguments. Let first C D, say v G C\D. Take f € Vy such that / = In <8> v - i.e., / = v: then / 6 C L 2tQ. H s 
and /| r = lr (8 v £ D L2( j-. H y Then 



C L2(n;H)X D L»<r;H) ~ U (g) P D y 
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i.e., (In <8> v)\r ^ lr <8> Pdv and accordingly Pc L2 , n . H) x.D L2 , r . H) f $ Vy. The case of D <£ C can be treated 
likewise. □ 

We mention that domination of semigroups can also be discussed. E.g., the following can be shown mimicking 
the proof of [59] Cor. 2.22]. This results extends [9j Prop. 2.8] and [55j Prop. 4.2]. 

Proposition 3.8. Impose the Assumptions [X3I and let Py be a positive operator. Let Si,S-z be L°°{T;C S (H))- 
function^. Define operators Si , 52 by 

Sig — Si o g and 6>2g = S2 ° g for all g G H^(T; H). 

Consider two sesquilinear forms a\,a2 defined by 

ai(f,g) := / (Vf(x)\Vg(x)) Hn dx + 1 f (V/(z)|V 5 (*)) H „_ a da(z) + (SMgirh 



and 

a 2 (f,g) := / (Vf(x)\Vg{x)) Hn dx + 1 f (V f(z)\Vg{z)) Hn . 1 da(z) + (S 2 f\r\9\r)y , 
Jn Jr 
both defined on Vy, and the associated operators Ay^j, Ay;s 2 . Then the following assertions hold. 

(1) The semigroup (e tAY - s i) t >o is dominated by (e* v,S2 )t>Oj *- e - 

\e tAY ^f(x,0\<e tA ^\f\(x,Z), t>0, /a 2 (l!xS;C), xeil^eZ, 
if and only if 

Re (S!fir\g\r)y > (S , 2|/||r||ff||r)y 
for all u, v G Vy such that uv > 0. 

(2) Let S\(z), S2(z) be positive operators for a.e. zeT. Then the semigroup (e tAY - s ^)t>o "is dominated by 
(e tAi, - s 2 ) t > if and only if S\(z) — S%{z) is a positive operator for a.e. zeT. 

Remark 3.9. In the usual theory of semigroup domination, both the dominating and the dominated semigroup 
have to act on the same space, or else one of them has to act on a space of scalar-valued functions, see |45) and 
references therein. This rules out several interesting case in our context, due to the fact the boundary conditions 
also determine the state space - and hence semigroups governing equations with different boundary conditions 
cannot been compared. E.g., it would be natural to expect that all semigroups (e tAy - s )t>o dominate the semigroup 
that governs the heat equation with (uncoupled) Dirichlet boundary conditions, provided that condition (|3.6I) holds. 

While it is known that many relevant properties are shared by the heat equation with either non-dynamic or 
dynamic boundary conditions, to the best of our knowledge a structural relation between these phenomena had 
not yet been observed. The following is a straightforward consequence of Theorem 13.61 and [T71 Prop. 4.3]. 



Corollary 3.10. Impose the Assumptions fX5l Then C£2 is left invariant under (e tA v s ) t>0 if and only if 
Cl 2 (q-,h) is left invariant under the semigroup governing the parabolic problem 

§- t u(t) = Au(t), t > 0, 

u(i) [r G Y, t > 0, 



(NDBC) 



^ + 5«(t)|r G Y±, t>0. 
u(0) = u , 



with time-independent boundary conditions. 

Observe that the semigroup governing (NDBC) is generated by the operator associated with 05 (with 7 = 0), 
but considered as a sesquilinear form acting on the Hilbert space {/ G H 1 ^; H) : /i r G y} ^ L 2 (VL; H) rather 
than Vy <-> £ 2 , cf. pi]. 



^ Here, we denote by L°°(T\ C S (H)) the space of all measurable and essentially bounded functions from T to C(H) with respect 
to the strong operator topology. 
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Example 3.11. As shown in [27, 9], remarkable properties of the (scalar-valued) heat equation with Wentzell- 
Robin (dynamic) boundary conditions include positivity and contractivity with respect to the co-norm of the 
semigroup that governs it. In the light of Corollarv \3.1(A these properties actually follow from the same properties 
enjoyed by the heat equation with corresponding Robin (time-independent) boundary conditions. 

Observe in particular that 

C p = C p y :=L p {fl;H) x (L P (T;H) fl y), p£ [l,oo], 

are Bochner spaces with respect to a suitable product measure. Assume both (e tAy ' s )t>o and its adjoint to be 
£°°-contractive: under the Assumptions 13.51 this can be characterized by means of Theorem 13. 6[ with Ch = 
(-oo, 1]h n [1, oo)h. 

Corollary 3.12. Assume both (e tAy - s )t>o and its adjoint to be C°° -contractive and let n > 2. Then (e tAy ' s )t>o 
extrapolates to a consistent family of operator semigroups on C p , p G [1, oo]. These semigroups are strongly 
continuous and analytic for p G (l,oo). 

Moreover, (e tAy - s )t>o is ultracontractive, i.e., it satisfies the estimate 

\\e tAy - s {\\ c ~ <M^t-t\\i\\ C 2 for all t€ (0,1], f G C 2 

where 

\[n- l,oo), if n > 3, 
" |(l,oo), if n = 2, 

and some constant M M . The same estimates are satisfied by the dual semigroup. 

Additional conditions ensuring strong continuity for p = 1 are known, cf. [6j §7.2.1] for the scalar case. 

Proof. The assertion on extrapolation follows applying a vector- valued version of Riesz-Thorin's interpolation 
theorem, cf. [331 p. 77]. The second assertion can be proved as in the scalar- valued case, applying a known 
characterisation of ultracontractivity (see [SJ §12.2]) based on standard Sobolev embeddings, cf. [53] Lemma 3.8]. 
It can be easily seen that all the involved techniques carry over to the vector-valued case. □ 

Remarks 3.13. 1) By [S2] Lemma 3.3], (e tAy - s )t>o consists of kernel operators for all t > 0. 

2) It is remarkable that the above mentioned criterion for ultracontractivity based on Sobolev embeddings 
only applies if n > 1 . In the scalar case, a common workaround is to deduce ultracontractivity from the Nash 
inequality. Unfortunately, the Nash inequality seems to extend to the vector-valued case only if the space H is 
finite dimensional. This is why we are not able to prove the above result in the case of n = 1 - which in particular 
corresponds to the relevant case of networks with infinitely many edges. 

A semigroup on an L 2 -space is said to be irreducible if the only closed ideals of L 2 left invariant under the 
semigroup are the trivial ones. If Y is a closed ideal of H, then clearly (e* Ay "°)t>o leaves invariant L 2 (fl;Y) x 
L 2 (T;Y), which is a closed ideal of C 2 . Thus, uncoupled boundary conditions jeopardize irreducibility. 

More generally, we observe that if V : fl — > C(H) is a strongly measurable function such that V(x) is an 
orthogonal projection onto a closed ideal of H for a.e. x € Q, then the subspace 

(3.7) I v := {/ G L 2 (fl; H) : f(x) G Range^x) for a.e. x G fl} 

is a closed ideal of L 2 (fl; H), too. In fact, all closed ideals of L 2 (fl; H) are of this form, as it is proven in |18) . 
Similarly, if the Assumptions 13.51 hold one can see that each closed ideals of C 2 is the range of an operator- valued 
strongly measurable mapping V defined on the product measure space fl © T and such that 

• V(x) is an orthogonal projection onto a closed ideal of H for a.e. x G fl and 

• Viz) is an orthogonal projection onto a closed ideal of Y for a.e. z G L. 

Proposition 3.14. Impose the Assumvtions HOI Then (e tAY -")t>o is irreducible if and only if Py is irreducible 
and fl is connected. 
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Observe that in the scalar case H = C the orthogonal projections on both subspaces of H are irreducible. 

Proof. It is clear that the semigroup is not irreducible if Q is unconnected, since it lets invariant the closed ideals 
consisting of those functions supported in any of the connected components. 

Let now Py be non-irreducible, i.e., consider a non-trivial closed ideal Jh of H such that PyJn C Jh- Then 
by Theorem l3.6l we conclude that Jl 2 {Vl\H) x is a closed ideal of L 2 that is left invariant under the semigroup, 
i.e., (e* Ai, <°) t >o is not irreducible. 

Let conversely (e* y '°)t>o be non-irreducible. Then there exists a non-trivial closed ideal of C 2 that is invariant 
under (e* Ar ' ) t >o. By Proposition 13.71 such an ideal is necessarily of the form Cl^(q.-h) x C^ct-h)- Now, we can 
apply Theorem 13.61 and deduce the claim. □ 

Remark 3.15. In the scalar case, it is known that irreducibility is equivalent to a strong parabolic maximum 
principle, provided that the semigroup is positive, cf. |59l §2.2] - but this characterisation fails to hold in the 
general vector-valued case. E.g., the heat semigroup (e tA )t>o on L 2 (K;R 2 ) is not irreducible because L 2 ($L;M. x 
{0}) is a non-trivial closed ideal left invariant under the semigroup. However, it does map nonzero positive 
functions f to functions e tA f satisfying e tA f(x) > for all t > and a.e. x G M. 

4. An example: Diffusion on a star-shaped network 

Throughout this section we consider the setting presented in Example 11.11 Observe that the Assumptions 13.51 
are satisfied whenever we discuss invariance of a set Ch which is either a subspace or an order interval containg 
0. We are going to present some interesting behaviour even in this elementary setting. Actually, same properties 
hold for more general diffusion on domains, rather than intervals. Also, by Corollary 13. 101 all the results in this 
section hold for the semigroups governing (NDBC) and (DBC) alike. Thus, we explicitly refer to the case of 
time- independent boundary conditions only. 

It has been proved in [T71 §5] that the semigroup governing (NDBC) is positive if Y = (1) (i.e., under so- 
called Kirchhoff boundary conditions) and not positive if Y = (1)^ (i.e., under so-called anti-Kirchhoff boundary 
conditions as considered e.g. in [33 |28j |6TJ [2] ) , provided that —S generates a positive semigroup on Y (i.e., —S 
is a real matrix with positive off-diagonal entries) . 

Similarly, assume that —S generates an i°°-contractive semigroup on Y and that H — C N . Then by [49l 
Lemma 6.1], this can be characterized by the fact that the entries s,y of S satisfy 

|sy| < Kesu for all i, 

cf. also [T71 Rem. 3.8.(2)]. Then one can prove that the heat semigroup is L°°-contractive under Kirchhoff 
boundary conditions for all feN, whereas in the anti-Kirchhoff case it is L°°-contractive if and only if N = 2. 

For the sake of simplicity, in the remainder of this section we let S = 0. 

A semi-infinite star with two edges can be identified with a line. More precisely, up to the canonical isometric 
isomorphism U defined by 

(Uf)(x) := ( f f ^ , x > 0, 

functions in L 2 (R; C) and in L 2 ((0, +oo); C 2 ) may be identified. Accordingly, a function (/i, f 2 ) £ L 2 ((0, +00); C 2 ) 
is called even (resp., odd) if fi(x) = f 2 {x) (resp., if fi{x) + f 2 {x) — 0) for a.e. x G (0, +00). More generally, we 
call a function / € L 2 (il;C N ) even (resp., odd) if f(x) € (1) (resp., if f(x) € (1 )) for a.e. x e ft. 

By Theorem l3.6l both the diffusion semigroups with Kirchhoff (i.e., Y = (1)) and anti-Kirchhoff (i.e., Y = (1)^) 
boundary conditions leave invariant the set of even functions as well as the set of odd ones. If N = 2, then it is 
easy to see that these are in fact the only boundary conditions leading to invariance of any of these both sets. 



I.e., e t ^f(x) is a nonzero, positive vector of R 2 . 
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Now, consider a semi-infinite star with only two edges, i.e., H = C 2 . Neglecting the trivial (uncoupled) 
boundary conditions defined by Y — {0} and Y = C 2 we can consider all 1-dimensional subspaces y = Y% of C 2 
by means of the parametrisation 

Py ( :=(. C ft Si ^ 2 CO f SC V ee[0,*), 
H \sm£ cos£ surf y s 

where denotes the range of the orthogonal projection Py ? . Observe that £ = 0, an( i £ = 

^ correspond to uncoupled Dirichlet /Neumann, to Kirchhoff, to uncoupled Neumann/Dirichlet and to anti- 
Kirchhoff boundary conditions, respectively, as can be checked directly. 

We are going to discuss the submarkovian property of the semigroup associated with these subspaces in 
dependence of £. A direct computation shows that the semigroup (e tAy « ,0 )t>o is positive if and only if £ € [0, §]. 
Furthermore, by Theorem 3.6 the semigroup that governs (NDBC) is L°°(il x S; C)-contractive if and only if 
Py t is L°°(S; (C)-contractive , i.e., if and only if the inequalities 

cos 2 £ + | sin£ cos£| < 1 and | sin£ cos£| + sin 2 £ < 1 

hold simultaneously. The former (resp., the latter) inequality holds if and only if £ ^ (0, j) U (^p,7r) (resp., if 
andonlyif^(f,f),(f,f )). 



















JT 


-JT 







cos 2 (£) + |sm(£) oos(£)| 
isin 2 (0 4- |s<»(£) ccs(e)| 



Figure 1 

Therefore, the L°°-contractivity of the semigroup associated with Kirchhoff boundary conditions represents a 
singularity. In particular, a submarkovian semigroup is generated exactly in the following five cases: 

• with uncoupled Dirichlct/Dirichlet boundary conditions, 

• with uncoupled Neumann/Neumann boundary conditions, 

• with uncoupled Dirichlet/Neumann boundary conditions, 

• with uncoupled Neumann/Dirichlet boundary conditions and finally 

• with Kirchhoff boundary conditions. 

Similarly, we can consider general boundary conditions denned by 1-dimensional subspaces of H for a semi- 
infinite star with 3 edges (H — C 3 ). They can be investigated by means of spherical boundary conditions, i.e., 
considering spaces y = Y^ j, that are ranges of the orthogonal projections 



(4.1) 



Py 



Py 



sin £ cos 2 <j> 
sin 2 £ sin <ft cos < 



sin £ sin <fi cos < 
sin 2 £ sin 2 <f> 



sin£ cos£ cos</> sin£ cos£ sin< 



sinf cos£ cosi 
sin£ cos£ sin< 
cos 2 £ 



£,0G [0,2tt) 



Analysing the behaviour of Py ( ^ in dependence of £, 4> as done above for Py i is less elementary. While the matrix 
is clearly positive if and only if £, 4> G [0, U [ir, it is not clear how to determine all the values £, </> leading 
to L°°-contractivity, i.e., all the values £, <j> such that the three functions 

2 



sin 2 £ cos 2 . 



| sin £ sin</> cos<f>\ + |sin£ cos£ cos0|, 

| sin 2 £ sin0 cos <f>\ + sin 2 £ sin 2 4-|sin£ cos£ sin</>| and £, <f> g [0, 7r) 

|sin£ cos£ cos0| + |sin£ cos£ sin</>| 4- cos 2 £ 

are simultaneously < 1, corresponding to the three conditions for i°°-contractivity associated with the three 
rows of the matrix Py ( . in (14.11) . 
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Figure 2B 



In Figure 2A we have plottecQ the level lines of the above functions for the value 1 (in violet, blue and red, 
respectively). This suggests that the ten parameter choices 

= (!,§), 

• 
• 



(£,4>) = (arctan v^f), 
(£, 4>) = - arctan \/2, f ), 
(£,0) = (arctan A^), 
(£» = (^-arctanv^,^). 



/ 37T 7T \ 

V 4 j 2 /> 

/ 7T 7T \ 

V 2 ' 4 /' 

/ 7T 37T \ 

V 2 ' 4 /> 

(5,0), 
(¥,o), 



lead to an L°°-contractive semigroup - as in fact can be checked directly. 
Observe that Y^ identifies Kirchhoff boundary conditions if and only if 

'i i r 



P^ 



i.e. , if and only if £ = arctan -\/2 and (f> = j . A direct computation shows that the remaining nine cases correspond 
to boundary conditions defined by means of spaces Y given by 

{(c, — c, — c) : c £ C}, {(c, c, — c) : c G 
{(c,c,0) : c£ C}, {(c,-c,0) : ce C 
{(0,c,-c) : ce C}, {(c,0,c):ceC 
While the last six subspaces only describe some decoupling of any of the three edges, we cannot find any physical 
interpretation for the first three boundary conditions. One can see that analogous boundary conditions give rise 
to L°°-contractive semigroups also in higher dimensional spaces H = C N for any A € N. 

It ought to be remarked that not all relevant values become evident through the above plot: one can see 
that decoupled boundary conditions arise with Yo,0 and Yj^ for all cf> £ [0,7r) as well as with Y^o and Y^z. for 
all £ £ [0,7r). Hence, using again the computations performed in the case of H = C 2 , we see that Yk^ lead 



{(c,~c,c):ceC} 
{(0,c,c):ceC}, 
{(c,0,-c) : ce C}. 



to i°°-contractivity for 4> e {0, f , § , and so do Y^ and Y^* for £ € {0, 4 



S [0,7r). We do not know whether further pairs (£,0) leading to L°°-contractivity exist. 



^ f , } as well as Y o>0 for all 



^ The figure has been obtained using Gnuplot 4.2 with a grid density of 1000 on both axes. For reference we have plotted the 
level lines for the value 0.5, too. On the £-axis (horizontal) we have highlighted the values ^, arctan \/2, ^ and n — arctan y/2. On 
the 0-axis (vertical) we have highlighted the values ^, ^ and 
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Moreover, a straightforward computation shows that Py ( + is a positive matrix if and only if (£, </>) G [0, -|] 
[0, Again, KirchhofT boundary conditions are a singularity in a "sea" of non-submarkovian behaviours. 
A similar procedure identifies all the 2-dimensional subspaces of C 3 , i.e., all ranges of the matrices Id — Py ( 



€,<#>> 



£, 4> € [0, 27r). However, plotting the level lines of the corresponding three functions (as we have done in Figure 
2B in violet, blue and red, respectively), does not suggest any new pairs (£, 4>) that lead to L°°-contractivity. 

The general case of a semi-infinite star with arbitrarily (finitely) many edges can be treated likewise, using 
known formulae for hyperspherical coordinates. 

As already remarked, the above results carry over to case of dynamic boundary conditions and should be 
compared with the known properties of the heat equation with Wentzell-Robin boundary conditions in the 
scalar case, cf. [lOl [55] and references therein. 

5. Dynamic boundary conditions on the normal derivative 

In this section we consider a different setting by discussing a new kind of dynamics on the boundary. While 
the dynamic boundary conditions introduced in (AV) involve the trace, dynamic boundary conditions on the 
normal derivative have also been considered in the literature, although less commonly (see |20ill6j ). Accordingly, 
the similar but different abstract initial-boundary value problem 



(AVN) 



■§- t u(t) = Au(t), t > 0, 

|^ g y, t>0, 

Su u (t)\r = SPyu ]r (t) + Py( 7 A r -S)^(t), t>0, 

u(0) = u , 

&(0)|r = «o, 



can be studied for 7 G R+ and S G C(L 2 (Sl; H)). The parameter S G C will be shown to influence the behaviour 
of the solutions to (AVN) in a curious way. 
Consider a sesquilinear form 65 defined by 

7i\ (91 



r i-\ i! — 1 (^fi(x)\^gi{x)) Hn dx-6(f 1 \ r \g 2 ) v da(z) 

HI \92j J J a 



{f2\9i\v)yda{z)+ 1 / {Vh{z)\Vg 2 {z)) H ^do-{z) + {Sf 2 \g 2 )y, 



with dense domain 



Wy := H 1 (SI; H) x (H s (T; H)ny), 
where s = 0if7 = 0ors = lif7>0. Mimicking the proof of Theorem l2.3l we deduce a corresponding generation 
result (cf. also the discussion in [THl §4.3]). 

Theorem 5.1. For any 7 € R + , S G C and S G C(L 2 (T; H)) the sesquilinear form by is continuous and elliptic 
(with respect to C 2 ). The operator By s associated with (65, Wy) generates an analytic semigroup (e tBy )t>o on 
C 2 with angle ^. The semigroup is compact if and only if Sl,T have finite measure, provided that H is finite 
dimensional. Moreover, by is accretive if 5 = — 1 and S is accretive; it is symmetric if and only if S = 1 and S 
is self-adjoint. In these cases the semigroups is contractive and self-adjoint, respectively. 



With a proof similar to that of Proposition 12.41 we can show the following, see also [15l §1.8]. 

Proposition 5.2. Assume SI to have C 2 -boundary. For any 7 G M+, S G C and S G C(L 2 (St; H)) the operator 
By.s associated with (bs,Wy) is given by 

•M r- W.. . A„, r- T2/r>. m A„._ r- T^lr. TT\ „„ J ^/ r- T 2 , 



D(By^ s ) = { ( J \ G Wy : Au G L 2 (Sl;H), Au r G L Z (T;H), and -± G L 2 (T;H) 



B 



y.s 



(A 

\5PyT Py^Av-S) 
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where T denotes the trace operator from to Hi (T; H), cf. [17, §7.1]. 

Thus, the semigroup associated with Bys yields the solution to the abstract initial-boundary value problem 

r fe(t) = Au(t), t>o, 
—ft) g y, t> o, 

< *£(t) = 5(P y u(t)) + (( 7 A r -<S)fK*)), t>0, 

w(0) = uo) 
I &(0) = 

Ouhabaz's criterion may be promptly applied to this setting, too. We omit the easy proof. 

Proposition 5.3. Impose Assumptions \3.5l Let C'h be a closed subspace or a closed order interval of H . 
Consider the closed convex subsets C_,2($-j ; __) an d Cy. Then (e tBy - s ) leaves invariant Cx2/Q ; m x Cy if and only 
if the compatibility condition 

SRc ( p c LHn ., a) h\r\ (I Pcy )h) y + Re (P cy h\ (i - Pc L2(ihH) ) h\r) y < Re (SP Cy f 2 \ (I P C ») h) y 

holds for all f\ £ i? 1 (f2; H) and all /_ £ y. 

Example 5.4. Impose Assumvtions \37h\ Then, by linearity (e tBY ' s ) is positive if and only if S — 1 and 

Re(SP Dy f 2 \(I-P Dy )f 2 ) y >0, 

i.e., if and only if 8 — 1 and the semigroup on H generated by —S is positive. 

Remark 5.5. It is easy to see that by similar methods one can also treat the parabolic problem 

jL u (t) = Au(i), t > 0, 

mPyu(t)\ dU = -Py™& + R lU (t)\ao, t>0, 

< &Py^ - -PyMt)\m + R*^, t>0, 
u(0) = u , 
w(0)|r = ui, 
V %r(0)|r = u 2 , 

for some i?i £ C(H x (Cl; H),y) and R 2 £ C{H l (n;H),y^). In this case the state space is L 2 (Vt;H) x}/x y 1 - . 
We omit the details. 
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